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We present a new method to solve the dynamics of disordered spin systems 
on finite time-scales. It involves a closed driven diffusion equation for the joint 
spin-field distribution, with time-dependent coefficients described by a dynamical 
replica theory which, in the case of detailed balance, incorporates equilibrium 
replica theory as a stationary state. The theory is exact in various limits. We 
apply our theory to both the symmetric- and the non-symmetric Sherrington- 
Kirkpatrick spin-glass, and show that it describes the (numerical) experiments 
very well. 
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Recently it has become clear JD that even mean-field models exhibit the ageing phenomena, 
familiar from experiments on real spin-glass that where hitherto assumed to be typical for 
short-range systems. This has led to a renewed interest in dynamical studies of mean-field 
spin-glass models and to valuable new insights into spin-glass dynamics away from equilibrium, 
see e.g. [||. In this letter we present a novel approach to analysing the dynamics of spin- 
glass models on finite time-scales, leading to a dynamical replica theory, which, in the case of 
detailed balance, incorporates equilibrium replica theory as a stationary state (including replica 
symmetry breaking, if it occurs). The formalism is built on a closure procedure which which 
we obtain a closed diffusion equation for the joint spin-field distribution. It constitutes the 
fixed-point of a series of previous dynamical studies . Our theory is proven to be exact for 
short times and in equilibrium. For intermediate times we can prove that it is exact if in the 
thermodynamic limit the joint spin-field distribution indeed obeys a closed dynamic equation. 
Here we discuss only the underlying physical ideas and the results of applying our theory to 



both the symmetric- and the non-symmetric 
mathematical details will be published elsewhere |£ | . 
and (simulation) experiment to be quite convincing. 



Sherrington-Kirkpatrick spin-glass. Full 
We believe the agreement between theory 



The generalised (asymmetric) version of the Sherrington-Kirkpatrick (SK) model 0, intro- 
duced in consists of N Ising spins Oi G {—1, 1} with infinite-range interactions Jy: 
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For i < j each of the random quantities Xij and , representing quenched disorder, are drawn 
independently from a Gaussian distribution with zero mean and unit variance. The evolution 
in time of the microscopic probability distribution Pt{&) is given by the master equation 
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in which Fk is a spin-flip operator Fk<fr(cr) = $(cti, . . . , — a^, . . . , ajsr) and the transition rates 
Wk(cr) and the local alignment fields hi(<r) are 



w k(s) = g I 1 ~ tanh[/3/i fc (<r)]] ^(ct) = J^cr, + (3) 

where (3 = 1/T is the inverse temperature. The mixing angle w € [0, 7r] controls the degree of 
symmetry of the interactions (|IJ). For a; = we recover the original SK spin-glass model 
Now the interactions are symmetric, the dynamics obeys detailed balance and (0) reduces to a 
Glauber dynamics, leading asymptotically to the Boltzmann equilibrium distribution Poo(c) ~ 
exp|— 0H(cr)] with the conventional Hamiltonian 

H(er) =-^T (TiJijCTj -0^2<?i (4) 

i<j i 

For to > 0, however, detailed balance is violated and equilibrium statistical mechanics no longer 
applies. For to — n the interaction matrix is fully anti-symmetric. 

For any given set of i macroscopic observables fi(cr) = (ili(er), . . . , O^(cr)) we can derive a 
macroscopic stochastic equation in the form of a Kramers-Moyal expansion. For determinis- 
tically evolving observables (in the thermodynamic limit) on finite time-scales only the first 
(Liouville) term in this expansion will survive, leading to the deterministic flow equation 

± n = v £ g pt(o-)*["-"(g)]Ei";i(g) [n(F ig )-n( g )] 
dt 1 n™oo 52<rPt(ir)s [n-n(*)] 

There are two natural ways for (||) to close. Firstly, by the argument of the subshell average in (||) 
depending on cr only through f2(er) (now Pt((r) will drop out), and secondly by the microscopic 
dynamics @ allowing for equipartitioning solutions (where p*(c) depends on a only through 
ri(er)). In both cases the correct flow equation are obtained upon simply eliminating pt(cr) from 
(||). We now make two assumptions: (i) the observables f2(er) are self-averaging with respect 
to the microscopic realisation of the disorder {xij, J/y}, at any time, and (ii) in evaluating 
the sub-shell average we assume equipartitioning of probability within the r2-subshell of the 
ensemble. As a result the macroscopic equation (^|) is replaced by a closed one, from which the 
unpleasant fraction is subsequently removed with a replica identity (see e.g. flQ] ): 

|nt = Jto o UmE-EE\ ,fl <( <Tl ) [n(V)-n((r l )]n*[n-n(0]) (6) 

"' ^ cr 1 cr™ i \ a =l I { x , y } 

For observables truly governed by a closed equation our closure procedure reduces to the natural 
one (in the sense discussed above), so we know by construction: if a closed self- averaging 
equation for flt(cr) exists, it must indeed be (^|). For the set of observables Jl(er) we now 
choose the (infinite dimensional) joint spin-field distribution: 

D(s,h;a) = ^^2S qt<rt S[h-hi(tr)] (7) 

i 

with the local fields (||). Since both the magnetisation m = jfY^i a i an d the energy (^) can 
be written as integrals over D(<j, /i;<r), the theory will automatically be exact (i) in the limit 
J — * 0, (ii) for t — > (upon choosing appropriate initial conditions) and (in) in the limit t — > oo 
(for systems evolving towards equilibrium) . To circumvent technical subtleties we assume that 
the distribution (Q) is sufficiently well behaved: we evaluate D t (s,h) in a number I of field 
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arguments and take the limit i — ► oo after the limit N — ► oo. A closed diffusion equation for 
D t (<;, ft) was also derived in JO]. Although similar in spirit to ours, the latter study employed 
a different closure procedure, lacking the properties of the present one of built-in exactness in 
various limits. 

We can now run the familiar machinery of replica theory and evaluate ([|) for the choice 
(0). The distribution Dt(s, ft) can be shown to indeed evolve deterministically. Details of this 
exercise, as usual involving a saddle-point problem, can be found in ||. The result is the 
diffusion equation 

h) = ~ [l+aanh(/3/i)] AK ft) - ~ [l-<rtanh(/3ft)] D t (?, ft) 

+ ^ | Afc,ft) [ft-0-J O <tanhC&ff)) Dt ] + + </ 2 [l-^tanh^i/))^] ^A(<T,ft) 

(8) 

with the short-hand (/(cr, H))r> = Yla I ^H D(<j, H)f(a, H). We find all interesting physics to 
be concentrated in a single driving term A: 

A[s, h;D]=- lim VOT 1 )^ { (tanh(/3 J ffi)cr a ) M (6[h-Hi]S^ ai [Hp -6-J Q m + iJ 2 cos wfflV)^) M 

a/3 

+ cosw (5[ft— ifi]5 C)<T1 cr a )jjf (tanh(/3i?i)[iJ / g — 6*— J m + iJ 2 coscj(J2^ er)^])^/ | (9) 
This involves an effective measure M in replica space: 

W H > a » M - JdffE ff M[fl lff ] 

M[Jf, cr] = exp < — irh • cr — -^J 2(T • Q a — i ^ D a (a a , H a ) 

- ^2 \H 6 J m+iJ 2 coswHV] ■ q~ l [H-d- J m+iJ 2 coswBVjj 



(10) 



with if = (Hi, . . . , i? n ) and cr = (cri, . . . , a n ). The n x n matrices {q, Q, R}, the n- vectors 
{m, m} and the functions D a (a, H) are obtained by extremisation of the surface 'J: 

m a m a + t;J 2 ^^ [QapQaj3 + CO8UlR a 0Rj3 a ] 
aa a a/3 

- i log det q + log [dH^ M [H,a] (11) 

2 J cr 

For the detailed balance case w = (the SK JtJ model) , the equilibrium state calculated within 
equilibrium statistical mechanics |L2j, is found to define a stationary state of (q). Here we 
restrict ourselves to replica-symmetric (RS) saddle-points of ^, an analysis of (nT) involving 
broken replica symmetry (RSB) a la Parisi |Q is the subject of a future study |13"| . It is 
a straightforward bookkeeping exercise to work out the RS saddle-point equations and the 
corresponding expression Ans for the driving term which, upon insertion into (^), controls 
the evolution of the joint spin-field distribution in RS approximation. Evaluating -4rs requires 
solving the RS saddle-point equations at each instance of time. In the usual manner one can also 
calculate the AT instability [1J] with respect to replicon fluctuations; here involving variation 
of three order parameter matrices, as opposed to one. 
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FIG. 1. Magnetisation m and energy per spin E in the SK model at T = 0, for J = (left) and 
Jo = 1 (right) . Solid lines: numerical simulations with N = 8000; dotted lines: result of solving the 
RS diffusion equation. 



J =0, im =0.3 



Jn=l, rXl n = Q.O 




FIG. 2. Field distribitions D t (cr, h) in the SK model at T = 0, for J = (left) and J = 1 (right) 
Histograms: numerical simulations with iV = 8000; lines: result of solving the RS diffusion equation. 
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FIG. 3. Magnetisation m and energy per spin E in the asymmetric SK model at T — 0, for Jo = 
(left) and J = 1 (right) . Solid lines: numerical simulations with N = 8000; dotted lines: result of 
solving the RS diffusion equation. 
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We test our theory by comparing the results of solving numerically equation (||) in RS ansatz, 
for the two model choices w = and uj — ^ir, with the results of performing numerical simula- 
tions of the discretised version of the stochastic dynamics (§,||) in a system of size N — 8000. 
Since solving (||) requires a significant computational effort, even within the RS ansatz, we 
restrict our experiments to zero external fields and to initial configurations with spin states 
chosen independently at random, given a required initial magnetisation. For the original SK 
model, obtained upon making the choice uj = 0, the results of confronting (||) with typical sim- 
ulation experiments at T = are shown in figures [l] and |2[ for Jq = (left pictures) and Jo = 1 
(right pictures). In figure [I] the top graphs represent the magnetisation m and the bottom 
graphs represent the energy per spin E\ for the two initial conditions too = and too = 0.3. 
Figure shows the corresponding distributions D t (a, h) for one particular choice of initial state 
(D t (l, h): upper graph in t = window, right graph in t > windows; D t (— 1, h): lower graph 
in t = window, left graph in t > windows). For Jo = 1 we were not able to calculate the 
solution of equation (g) up to t = 6, due to the critical behaviour of the saddle-point equations. 
For the fully asymmetric model, corresponding to uj = ^ir, one finds much simpler equations, 
due to a decoupling of the spins from the fields. In this case equation (||) in fact allows for 
distributions D t (c,, h) which remain of a Gaussian form at all times, in accordance with 15 



The results of confronting (^) with T = simulation experiments for the asymmetric uj = -^n 
SK model are shown in figures || and [|, for Jo = (left pictures) and Jo = 1 (right pictures); 
for the three initial conditions too = 0, mo = 0.3 and to = 0.6. Figure ^ shows the distributions 
D t (a, h) for one particular initial state (D t (l, h): upper graph in t = window, right graph in 
t > windows; D t (—1, h): lower graph in t — window, left graph in t > windows). 

In this letter be have shown how one can solve the dynamics of disordered spin systems 
on finite time-scales with a dynamical generalisation of equilibrium replica theory. Although 
we have restricted our analysis by making the replica-symmetric (RS) ansatz, on the time- 
scales considered the agreement between theory and simulation experiment is already quite 
satisfactory. At this stage we need more efficient numerical procedures in order to extend the 
time-scales for which we can solve the equations of the theory. This would enable us to compare, 
for instance, with data such as the ones in |Iq] , and to investigate the possible existence of 
stationary states other than the one corresponding to thermal equilibrium. Our theory is by 
construction exact in various limits. Its full exactness depends crucually on whether the joint 
spin-field distribution indeed obeys a closed self-averaging dynamic equation, which is difficult 
to verify. We plan to investigate several approaches to this problem in the near future. Firstly, 
we want to apply our formalism to disordered spin systems for which the dynamics has been 
solved by other means, like the toy model |J, or the spherical spin-glass [|l7|. Secondly we want 
to try to derive a diffusion equation for the joint spin-field distribution, starting from the exact 
equations for correlation- and responsefunctions, as obtained from the path-integral formalism. 
The latter approach involves (rather complicated) closed equations for the correlation function 
C(t,t') and the response function i?(i,i'), with two real-valued arguments each (two times). 
The present formalism also involves two functions D t (l,h) and D t (— 1, h), with two real- valued 
arguments each (one time and one field). It is therefore quite imaginable that both formalisms 
constitute exact discriptions of the dynamics of disordered spin models. 
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